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^^ ■ Abstract 



k>( \ We study some particular solutions to the Navier-Stokes-Poisson equations with density- 

C^ ' dependent viscosity and with pressure, in radial symmetry. With extension of the previous 

known blowup solutions for the Euler-Poisson equations / pressureless Navier-Stokes-Poisson 
with density-dependent viscosity, we constructed the corresponding analytical blowup solu- 
tions for the Navier-Stokes-Poisson Equations with density-dependent viscosity and with pres- 
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1 Introduction 

The evolution of a self-gravitating fluid can be formulated by the isentropic Euler-Poisson equations 
of the following form: 

Pt+V • (pu) =0, 
p\u^ + (u- V)u)] + VP = -(5pV$ + vis(p, u), (1) 

A$(t,a;) =a(A^)p- A, 

where a(N) is a constant related to the unit ball in R^ : a(l) = 2; a(2) = 27r and for A^ > 3, 

a(iV) = JV(JV - 2)T/o/(JV) ^ JV(JV - 2) ^^^^^ ^ ^^ , (2) 

where yoZ(A^) is the volume of the unit ball in R^ and F is a Gamma function. And as usual, 
p = p(i, x) and w = M(t, x) e R^ are the density and the velocity respectively. P = P{p) is the 
pressure. And A is the background constant. 

When (5 = 1, the system can model fluids that are self-gravitating , such as gaseous stars. In 
addition, the evolution of the simple cosmology can be modelled by the dust distribution with- 
out pressure term. This describes the stellar systems of collisionless and gravitational n-body 
systems j8] . And the pressureless Euler-Poisson equations can be derived from the Vlasov-Poisson- 
Boltzmann model with the zero mean free path fTU]. For A^ = 3, the equations ([T]) are the classical 
(non-relativistic) descriptions of a galaxy in astrophysics. See [2] and [3], for details about the 
systems. 

When 5 = — 1, the system is the compressible Euler-Poisson equations with repulsive forces. The 
equation ([ills is the Poisson equation through which the potential with repulsive forces is deter- 
mined by the density distribution of the electrons. In this case, the system can be viewed as a 
semiconductor model. See [5] and [6] for detailed analysis of the system. 

When (5 = 0, the potential forces are ignored. The system is called Euler / Navier-Stokes equations. 
See [1] , [12] and [T7] for detailed analysis of the system 

In the above system, the self-gravitational potential field $ = $(t, a;) is determined by the 
density p through the Poisson equation. 
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And vis{p,u) is the viscosity function: 

vis{p, u) := v(m(p) V -u). (3) 

In this article, we seek the radial solutions 

p(i, x) = p{t, r) and u = -V{t, r) =: -V , (4) 

r r 

1 /? 
with r = I 'Ylii=i ^i I ■ By the standard computation, the Euler-Poisson equations in radial 

symmetry can be written in the following form: 

N -1 

Pt + Vpr + pVr + pV = 0, 

r 

P iVt + VVr) + Pr{p) + 5p^r (p) - Mr (p) {^V + Vr) + p{p) [~ (^) V + ^F, + Kr] ■ 

(5) 
Under a common assumption for, the viscosity function can be defined: 

pip) := Kp^ (6) 

where k and 9 > are the constants. For the study of the above system, the readers may refer [T3], 
|16j . [18] and [21]. In particular, when 9 = 0, that returns the expression for the only V^-dependent 

viscosity function: 

N- 1 

Pt + Vpr + pVr + pV = 0, 

(7) 

p {Vt + VVr) + Pr{p) + Sp'^r {p) = At [Kr + ^K- " (^) V] , 

which are the usual form of the Navier-Stokes-Poisson equations. The equations ([T])i and ^2 
(vis{p, u) ^ 0) are the compressible Navier-Stokes equations with forcing term. The equation 
(IT])3 is the Poisson equation through which the gravitational potential is determined by the density 
distribution of the density itself. Thus, the system ([T|) is called the Navier-Stokes-Poisson equations. 

Here, if the vis{p,u) = 0, the system is called the Euler-Poisson equations. In this case, the 
equations can be viewed as a prefect gas model. For iV = 3, ([1]) is a classical (nonrelativistic) 
description of a galaxy, in astrophysics. See [5], [TS] for details about the system. 

P = P{p) is the pressure. The 7-law can be applied on the pressure P{p), i.e. 

P{p) = Kp'':=^, (8) 
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which is a common hypothesis. The constant 7 = cp/c^ > 1, where cp, Cv are the specific heats 
per unit mass under constant pressure and constant volume respectively, is the ratio of the specific 
heats, that is, the adiabatic exponent in ([8]). With X = 0, we call that the system is pressureless. 

For the local existence of the Euler-Poisson equations, the interested reader may refer the results 
of Makino [M] , Gambin [9] and Bezard [1] . 

In the following, we always seek solutions in radial symmetry. Thus, the Poisson equation ([Ijs 
is transformed to 

r^-i$„ (t, x) + iN- 1) r^~Hr=a (N) pr^-\ (9) 



^r^^^ / {p{t.s)-k)s''-Hs. (10) 



In this paper, we concern about blowup solutions for the A^-dimensional pressureless Navier-Stokes- 
Poisson equations with the density-dependent viscosity. And our aim is to construct a family of 
blowup solutions. 

Historically in astrophysics, Goldreich and Weber [TT] constructed the analytical blowup (col- 
lapsing) solutions of the 3-dimensional Euler-Poisson equations for 7 = 4/3 and for the non-rotating 
gas spheres. After that, Makino [H] obtained the rigorously mathematical proof of the existence 
of such kind of blowup solutions. And in [7], Deng, Xiang and Yang extended the above blowup 
solutions in R^ (N > 3). Recently, Yuen obtained the blowup solutions in R^ with 7 = 1 in |.21j . 
The family of the analytical solutions are rewritten as: 
for A^ > 3 and 7 = (2A - 2)/N, in [7]: 

. , . «7^2/(^)"/("-^\ for r < am,, ^^^ ^ „(,) 



a{t) 
0, for a{t)Z^, <r. 

q(^)= ^/^xjv-i ' "('^) = Oo > 0' o.{Q)= ai, 
where p = [N{N - 2)A]/(2A - 2)K and the finite Z^ is the first zero of y{z); 



(11) 
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for TV = 2 and 7 = 1, in [gT] 



P{t.r) 



1 eV(r/a(t))^y^^^^-^^^ 



a{t) 



a(i)=— — ■, a(0) = flp > 0, a(0)= oi, 
a[t) 

1 2tt y^^^ 

y{z) + -y{z)+^e = fi, y(0) = a, ?;(0) = 0, 
z K 



(12) 



where K > 0, ^ = 2X/K with a sufficiently small A and a are constants. 

For the construction of analytical solutions to the Euler, Navier-Stokes and pressureless Navier- 
Stokes-Poisson equations, readers may refer to the recent results in [3D], [TH], [23] and [31] . 

In this article, it is natural to extend the more general results to cover the full system ((SJ — - 
Navier-Stokes-Poisson equations with density-dependent viscosity and with pressure. In short, we 
successfully deduce the system ([5]) into the much simpler ordinary differential equations. 

Theorem 1 For the Navier-Stokes-Poisson equations in R'^ {N > 2) with 7 = ^^^ , ^ there 
exists a family of solutions, those are: 



p{t,r) = 



0, for a{t)Z^ < r. 



V{t,r)^^r 



(13) 



where f{z) and a{t) are the following functions and the finite Zf^ is the first zero of f{z): 

(1) with A = 0, and 

(la)9^^^, 

a{t) — mt -\- n, 
^Kf{zV-^f{z) - mK9Nfiz)0-^f{z) + ^ J^ f{s)s^-'ds, /(O) = a > 0, 

where m, n > and a are constant; 

(ih)e = ^^, 

a{t)^{mt + nf'^, 
jKf{zr-^f{z) - 2mK0fizy-^fiz) + ^ /; fis)s^-^ds = ^<^-^^"'% , /(O) = a > 0. 



(14) 



(15) 



In particular, for m < 0, the solutions |j^[ j and US]) blow up at the finite time T — ~m/r 
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(2) with 6A> and 6 = 2i^ , we have 



/Sc,(N)A, 



lit) 
(Kj ~ NSaiN)A) f{zY-^f{z) + ^ /„^ f{s)s^-Hs, /(O) = a, /(O) = a > 0. 



(16) 



2 Separable Solutions 



Before we present the proof of Theorem [1] Lemmas 3 and 12 of [35] are quoted here. 



Lemma 2 (Lemma 3 of |22p For the equation of conservation of mass in radial symmetry: 



N -1 

Pt + Vpr + pVr + pV = 0, 

r 



(17) 



there exist solutions, 



pU,.f±Mm V(tr)-^-^r 



(18) 



with the form f > e C^ and a{t) > G C^. 



It is clear to check our solutions to satisfy the Navier-Stokes-Poisson equations: 
Proof. As we use the functional structure of the above lemma, our solutions (|13p fit well to the 
mass equation fUi. 
(la) For the momentum equation ([5]) 2, we have: 



p{Vt + V-Vr)+K ip'>)^ - K{p'')r [ V + Vr)~ Kpr{Vrr + K —V) (19) 



^ Jo 



P^^ + K 
a{t) 



J\a(t)) 






(20) 






ait) ^ r^-^ Jo a{t) 



-fK 



■film) 



7-1 



d_ 
dr 



/(^ 



W)' 



a{t)^ 
6a{N) /"■/(^) ^_i 



a{t)^ 



kO 



a{t)^ 



d_ 
dr 



■lya{t}' 

^{tj^ 



Na{t) 
<t) 



r^-i Jo a{t) 



\N 



ds. 



(22) 
(23) 
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with a{ 


t) = mt + n. 










Then, we have: 








= jK 




7-1 




-k9 







Nd{t) + p- 



5a(N) r f{^) 



.Af-l 



,N-1 



lit) 



N 



ds 



(24) 



'I a(t)^(7-2)+W+l a(i)iV(e-2)+W+2 



^W + 7^X-^^""^'1 ^''^ 



P jiKfi^V-'fi^) ^oNfi^y-^fi^) sa{N) r/^^'^ m ^_, 



a(i)JV-l I (j(i)Ar(7-2)+2 (j(i)Ar(e-2)+3 



d(t) 



a{t) 






a{tr 



s"~'ds' 



(26) 



And we plug with the condition for 7 = ^^q^ and 9 = ^^ "^ and with the new variable z := r/a{t): 



-^ l^^KJ{zy-^f{z) - mneNf{z)<^-^j{z) + ^M /^ /(.).^-id. ^ . (27) 



In the theorem, we require the ordinary differential equation for y{z) is 



-/Kfizr-^fiz) - m^eNf{zf-\f{z) + ^ /; f{s)s^~'ds = 0, 



y{0) ^a>0. 



Therefore, our solutions satisfy the momentum equation. 

(la) Similarly, for A = 0, 7 = ^^ and 6 = 2^, we have: 



(28) 



'A^-1 \ A^-1 A^-1 

p{Vt +V-Vr)+K {p^)^ - K{p'')r ( V + Vr)- Kpr{Vrr + V,. —V) (29) 



Sa{N) 

^ •.'0 



ps^'-^ds 



(30) 



a{t) 



jya(t)> 



a{t)^ 



Ndit) ^^SaiN) r f(M,.^^,, (31) 



ait) 



a(^) 



2.(2. -1)^2 
p , , r^— r + 7Ar 



Ke* 



(mt + n)^ 
a{t)^ 



/(^) 



a(t)^ 



7-1 



d_ 
dr 



a(i) 



N 



d_ 
dr 






2m{mt + n) + p- 



_i , Sa{N) r.f{^) jv-i 



r^-i ./o a{t)^ 



s''-'ds, 



(32) 



(33) 
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with a{t) = {mt + nf/'^. 
Then, wc have: 



2(2- N)m^r 

/9— ^ \~"iK 



N'^{mt + ny 



fia 



(t)' 



ait) 



N 



7-1 



ii-a 



it)' 



a(t) 



N+l 



2mK9 






^^"(*)^ imt + n)('-^y^+p^-4^ rl^s^-'ds 



ait) 



N+2 ' 



r~-i Jo ait) 



\N 



N'-'imt+n)'-' "' a(t)«(^-2) + w+i 



a(i)N(e-2) + iv+2 (,"'■'' T" "-^ 'T- ^jv-1 Jo a(t)N * "* 



2(2-JV)m'' (mt+«)^<"-i)/" , 7J^/(^)^'V(^) 



N'-' (rnt+n) 



^ r 



imt + n)^(^-^y^ 1 2mnef{^)0-y(^) 

a(t)N(e-2)+3 



{mt + n)^-^ + ^^^jyi 



a(t)"W-2)+2 



(ty 



ds 



2{2-NW r , 7'^/(^)^''/(^) 



ait) 



W-1 



(mt+n)»("(»-2)+3) 



(?7it + n) 



,-1 <5a(JV) rr/a(t) f (s) „Ar-lj„ 



;rwT^r^ 



a(t)' 



«(i) 



N-l 






(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



2Ar-2 



and 9 = 2^ for iV > 2. 



with 7 - AT — - - 2W 



and with the new variable z :~ r/ait): 



p r 2(2 - N)m-_^ ^ ,Km^-^m - 2m.em^-^m + ^ / /(.).--ci. 



In the theorem, we require the ordinary differential equation for /(z 



(40) 



SajN) r- f(a\aN-lj, _ 2{N-2)m^ 



jKf{z)''-\fiz) - 2mnef{z)<^-yiz) + '-^ J^ f{s)s^-'ds = - ^V ' z, 

2/(0) = a > 0. 



(41) 



Therefore, our solutions satisfy the momentum equation. 
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(2) for (5A > and 7 = 6* = ^^, we have 



'iV-1 \ iV-1 A^-l 

p{Vt + y • K) + if (p^), - k(/), ( V + Vr]- Kpr{Vrr + K- ^V) 






5a{N)K 
P^j^r + K 



fi-a 



W)' 



ait) 



N 



a{t)^ 



Nd{t) 



(42) 
(43) 

(44) 



<5a(7V) r/(^) jv-ij M^)A 



^iV-l 



z(t) 



Af 



s''-'ds~p- 



rN-1 



„JV-1 



ds. 



(45) 



/ 5a{N)A 

With a(t) = eV « , then we get: 



Kj 






7-1 



a(i)^+i 



k9 



a{tr 



e-1 



llf^./NS^^iNjA + p^"^^^ r^^^wKN-i 



a{t) 



r.N-1 



a{t) 



N 



ds 



(46) 



p<jK 



a(t)^ 



7-2 



J\a(t)' 

a{t)^+^ 



K9^yNSa{N)A 






/(^) , fo(iV) f /(^)^^„i^^ 



a(i)^+i r^-i 7o aW 



\N 



(47) 



f 7i^/(^)^-V(^) ^VA^fo(JV)A/(^)^-V(^) 
^1 a(t)^(7-2)+w+i a(i)iV(e-2)+iv+i 



fc(iv) r -^(^1,^-1^,^ 



r^-i 7o «(t) 



W 



(48) 



p f7^/(^r-V(^) «V^'5a(iV)A/(^)'^-V(^) , fo(7V) r/(^) „_, 



a(t)^-i j a(t)^(^-2)+2 



a(t)^(»-2)+2 



a(*) 



iv^n- "'0 



a(t) 



N 



■s"-'ds 



(49) 



P f7^/(^r-V(^) n9^N5a{N)Afi^r-'fi^) Sa^N) r'<'^ /(.) ^_, 



a{t)^-^ \ a(t)^(^-2)+ 



a(t)^( 



7-2)+2 



+ 



a{t) 



N-l 

■n^tt Jo 



^ ^"■«^>^"^<^> - ''^v^^^<^«^'""^(55' " ^ 



a(t) 



■N^T- JO 



o(t)^ 
(50) 

r/a{t) j^^ 



s"-'ds 



a{t)^ 



s^-irfs 



(51) 
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2N-2 



with 7 = 61 JJ-. 

We require the corresponding ordinary differential equations f{z) with z := r/a{t) 



^K - Ke^N5a{N)K) f{zr<-^f{z) + ^ j^ f{s)s^-Hs = 
y(0) = a > 0. 



(52) 



The proof is completed. ■ 

Remark 3 If the complex number solutions {p, u) e (7^+^ is considered, we have the corresponding 
solutions for 5A < 0: 



/(^) 



Pit,r)^^^^, V{t,r)^t^^r, 

ait) = e^V^\ (53) 

(7if - iK9^-NSa{N)A) f{zy-^f(z) + ^ /^^ f{s)s^-'ds = 0, /(O) - a, 
where i is the complex constant. 

Remark 4 Our method can be easily extended to the corresponding systems with linear damping: 



Pt+V ■ [pu) =0, 
p[mj + {u- V)u] + Ppu + VP = -5pV^ + vis{p,u), 



(54) 



A$(t,a;) ^a{N)p- K, 



where P > is a constant. 
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